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Observational signatures of sub-photospheric radiation mediated 
shocks in the prompt phase of GRBs 
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A shock that form below the photosphere of a GRB outflow is mediated by 
Compton scattering of radiation advected into the shock by the upstream fluid. 
The characteristic scale of such a shock, a few Thomson depths, is larger than any 
kinetic scale involved by several orders of magnitudes, hence, unlike collisionless 
shocks, radiation mediated shocks cannot accelerate particles to non-thermal 
energies. The spectrum emitted by a shock that emerges from the photosphere 
of a GRB jet, reflects the temperature profile downstream of the shock, with 
a possible contribution at the highest energies from the shock transition layer 
itself. We study the properties of radiation mediated shocks that form during 
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the prompt phase of GRBs, and compute the time integrated spectrum emitted 

by the shocked fluid following shock breakout. We show that the time integrated 

emission from a single shock exhibits a prominent thermal peak, with the location 

|>- ■ of the peak depending on the shock velocity profile. We also point out that 
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1. Introduction 

There is mounting evidence that a significant fraction of the bulk energy of GRB outflows 
is dissipated below the photosphere (e.g., Eichler & Levinson 2000; Peer, et al. 2006, Ryde 
& Peer 2009; Peer, et al. 2012). Such dissipation may be accomplished through formation 
of shocks in a hydrodynamic flow (Levinson & Eichler 1993; Eichler 1994; Meszaros & Rees 
2000; Giannios 2011), magnetic reconnection in a magnetically dominated flow (Lyutikov 
& Blandford 2003; Giannios & Spruit 2007; McKinney & Uzdensky 2012), and perhaps 
other means (e.g., Beloborodov 2010). Recent simulations of jet breakout from the stellar 
envelope of a progenitor star indicate that strong tangential collimation shocks form below 



multiple shock emission can produce a spectrum that mimics a Band spectrum. 
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the photosphere of the expanding outflow, in a region of moderate optical depth (1 < t < 50) 
(Lazzati et al. 2009). The photosphere is located at a distance of several hundreds stellar 
radii from the central engine, r p h ~ 10 12 — 10 13 cm, with the exact location depending 
on the type of the progenitor star. These collimation shocks can dissipate a considerable 
fraction of the bulk energy, leading to variable emission on intermediate timescales (Morsony 
et al. 2010). In addition to the collimation shocks, intermittencies of the central engine are 
expected to lead to formation of internal shocks, that are commonly thought to be the origin 
of the short-timescale variability (down to milliseconds) observed in many bursts. These 
shocks can be produced above the photosphere or below it, depending on the bulk Lorentz 
factor of the expelled outflow and the duty cycle of the central engine (Bromberg et al. 2011, 
BML11). It has been shown (BML11) that for moderate bulk Lorentz factors, T < 500, and 
a reasonable duty cycle these shocks should form in a region of moderate optical depth, 
1 < r < 10 3 . 

Shocks that traverse an optically thick plasma are mediated by radiation when the shock 
velocity (measured henceforth in units of c) satisfies 

f3 s » 4 x lO-yf, (1) 

where n = 10 15 nis cm -3 is the density of the unshocked plasma (Weaver 1976). This condi- 
tion is always satisfied for shocks that form below the photosphere of a GRB outflow during 
the prompt phase, hence they must be mediated by radiation. In case of a nonrelativistic ra- 
diation mediated shock (NRRMS), f3 s << 1, the length scale of the shock transition layer 
can be estimated by equating the photon diffusion time across the shock, t' D ~ n e OT^ /c, 
with the shock crossing time, t' s = A' s /(cj3 s ). This yields A' s ~ (na T fis)~ l - In case of a 
relativistic radiation mediated shock (RRMS) that propagates in a cold medium, the shock 
width, as measured in the shock frame, is roughly ~ 0.1(n<7r) -1 7 2 for sufficiently large 7 U 
(Budnik et al. 2010 (BKSW10); Nakar and Sari 2012). Here 7„ is the Lorentz factor of the 
upstream fluid with respect to the shock frame. The width of a RRMS which is dominated 
by photon advection, as discussed below, may have a different scaling, but is also on the 
order of (o^n) -1 ( Levinson and Bromberg 2008). Thus, the width of radiation mediated 
shocks span a range between a few to a few tens Thomson depths. 

The structure and emission of a RRMS were computed in BKSW10, assuming that the 
upstream plasma is cold. Under this condition the dominant photon source inside and just 
downstream of the shock transition layer is Bremsstrahlung emission by hot electrons and 
positrons that are produced via annihilation of counter streaming photons. Since, as will be 
shown below, the photon generation rate is too low to allow for the radiation produced inside 
the shock to thermalize, the characteristic temperature in the immediate post shock region 
largely exceeds the black body limit. Typically, 60 keV ^ T s <, 200 keV for f3 s > 0.6 (Katz 
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et al. 2009 (KBW09)). The spectrum of the radiation produced inside the shock exhibits a 
tail that extends to energy in excess of , y u m e c 2 in the shock frame (owing to pair production 
in the near upstream region). For typical GRB parameters (as well as for other systems) 
the length scale of the shock transition layer, A' s ~ 10 9 n^ cm, is larger by many orders 
of magnitudes than any kinetic scale involved, and in particular the plasma skin depth, 
c/co p ~ l^i5^ 2 cm > an d the gyroradius of protons, tl ~ 3^/3 (B/10 6 G) _1 cm. An immediate 
consequence is that radiation mediated shocks cannot accelerate particles, contrary to what is 
commonly assumed. Another consequence is that the hot electron-positron pairs are tightly 
coupled to the protons through electromagnetic interactions and, therefore, thermalize on 
scales much shorter than the Thomson mean free path. The local distribution of the electrons 
and positrons can therefore be taken to be Maxwellian everywhere. 

The conditions invoked in BKSW10 are anticipated in shocks that breakout of the 
envelope of an exploding star (e.g., Weaver 1976; Chevalier 1992; Rabinak & Waxman 2011; 
Nakar & Sari 2010, 2012; Sapir et al. 2011), and perhaps in some other situations. However, 
as shown in BML11, under the conditions prevail below the photosphere of a GRB outflow, 
the fluid upstream of a sub-photospheric shock advect photons at a high enough rate to render 
photon production in the shock transition layer negligible. As a result, the temperature in 
the immediate downstream is significantly reduced, but still exceeds the black body limit. 
The spectrum emitted by the hot plasma downstream of a shock that breaks out of the 
photosphere depends on the temperature profile behind the shock. In what follows, we 
analyze the properties of sub-photospheric shocks that form in the prompt phase of GRBs, 
and compute the spectrum emitted by a single shock that emerges from the photosphere. 
In §2 we examine the conditions behind a planary shock, both in the relativistic and in the 
nonrelativistic regimes. In §3 we apply the results obtained in §2 to analyze the properties of 
a radiation mediated shock propagating in a GRB jet, and to compute the emitted spectrum. 
We conclude in §4. 



2. Structure of a radiation mediated shock 

Consider a planar shock propagating in a homogeneous medium consisting of protons, 
electrons and possibly seed radiation. Under the conditions envisioned here, the shock is 
mediated by Compton scattering of radiation produced at about one Thomson depth down- 
stream of the shock transition layer - the deceleration zone where the flow velocity changes 
from its far upstream value to its terminal value downstream. In the frame of the shock, the 
jump conditions read: 
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WullPl +Pu = w d -f 2 d /3 2 d + p d , 



(3) 
(4) 



where subscripts u and d refer to the upstream and downstream values of the fluid parame- 
ters, respectively, rib denotes the baryon density, w the specific enthalpy, f3 the fluid velocity 
with respect to the shock frame, and 7 = (1 — /3 2 ) _1//2 the Lorentz factor. For the situations 
considered in this paper the pressure, when important, is always dominated by radiation. 
Thus, the specific enthalpy, both upstream and downstream of the shock transition layer, 
can be approximated as w = ribm p c 2 + 4p r /3, where p r denotes the radiation pressure. 



We consider first the ultra-relativistic limit of Eqs. (J2J)- (J4j) . For shocks that form well 
above the coasting radius of a GRB jet the ratio of radiation pressure and rest mass energy 
density upstream of the shock, p = p r u/( n bu m P c 2 ), is always small, as implied by Eq. 
below. Thus, to a good approximation we can neglect the pressure of the upstream fluid; 
that is, we set p u = in Eqs. ()3]) and (TjJ. Then, the jump conditions, augmented by the 
equation of state Wd = n^m^c 2 + 4p r ^/3 downstream, reduce to 



Now, since the upstream flow is supersonic, the photon density there, n rU} is always smaller 
than that required to establish a full thermodynamic equilibrium downstream. Therefore, 
photons must be produced inside and downstream of the shock transition layer in order 
for the radiation to reach thermal equilibrium. For the fast shocks considered here the 
thermalization time is much longer than the shock crossing time (see below). Consequently, 
the temperature T s in the immediate post shock region, where the shock velocity approaches 
its terminal value (3d, is well in excess of the black body temperature Td- Henceforth, we 
denote quantities associated with the radiation field in the immediate post shock region by 
subscript s, and quantities in the far downstream region by subscript d. The energy density of 
the radiation in the immediate post shock region, e rs , is then given by Eq. fl6]). For the infinite 
planar shock studied here the energy downstream is conserved, so that e r d = e rs . However, 
under the conditions considered in the following section, where the shock propagates in an 
inhomogenous, relativistically expanding flow, the radiation cools adiabatically during the 
evolution of the shock and e r d decreases with time. 



2.1. 



RRMS 



f3 d = 1/3, 
e rd = 2n bu jl(3lm p c 2 . 



(5) 
(6) 



The thermalization length depends on the rate at which photons are generated. For the 
situations addressed here the relevant photon sources downstream are free-free emission and 
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double Compton (DC) scattering (Svensson 1984) J^J The corresponding photon production 
rates per unit volume are denoted by riff and h DC , respectively, and are given explicitly 
in appendix [X] Sufficiently far downstream the radiation eventually thermalizes. Then, the 
photon density reaches the black body limit, and is given by n r d = e T d/kTd with e rd = aT^. 
Using Eq. (jSJ) one finds 

n rd = 2.8 x 10 21 n^ 5 ( 7t A) 3 / 2 cm" 3 , (7) 
kT^OAnl^M 1 / 2 keV, (8) 

where n bu = 10 15 n^ 15 cm -3 . To get a rough estimate of the thermalization length, we 
consider each process separately. When double Compton dominates, the thermalization 
time is t' DC = n r d/fiE>c- The thermalization length can be expressed as L' DC = (3dCt' DC . 
When measured in units of the Thomson length it reads 

T ' DC = a T n bd L' DC ~ 10 6 A^n^7-\ (9) 

where Eq. (IA3I) has been employed with A^c — lOA^ci an d T = Likewise, for thermal- 
ization by free- free emission we obtain from Eq. (IA2I) with Aff = lOAffi and Eq. (j7]), 

r' ff = a T n bd L' ff = 1.5 x lO^n^W 4 . (10) 

These rough estimates are in good agreement with the detailed calculations of the temper- 
ature profile outlined in appendix [S] (see figured]). Evidently, the length scale over which 
thermodynamic equilibrium is established largely exceeds the shock width. 

The temperature in the immediate post shock region, at a few Thomson depths down- 
stream of the shock transition, depends on the seed photon density. In what follows we 
focus on advection dominated shocks for which photon generation inside the shock (but 
not in the thermalization layer) can be neglected. Then, the photon density just down- 
stream of the shock transition layer is the advected density compressed by the shock: 
n rs = fiulun r u/{fidla) — V^lufiunru- Thermal Comptonization of the advected photons 
by the hot electrons in the immediate downstream results in a Wien spectrum with an av- 
erage photon energy < hv >= 3kT s = e rs /n rs = e rs / (■\/&'y u f3 u n ru ) (BKSW10, BML11). By 
employing fl5]) and we obtain 

kT s ~ 23n 4 "HA keV, (11) 



1 Synchrotron emission is another potential source of photons. However, the characteristic frequency of 
cyclotron photons emitted by thermal electrons is hv syn — 0.03£?6 eV, where Bq = Bj (10 6 G) ~ 1 is roughly 
the equipartition magnetic field strength of a typical GRB near the photosphere. It can be readily shown 
that such soft photons are self-absorbed well before their energy is increased significantly, and therefore 
cannot contribute to the thermal peak. 
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here n = n ru /rib u = 10 4 n4 is roughly the dimensionless entropy per baryon carried by the 
upstream flow. 

In the thermalization layer the temperature changes from T s to T d < T s , with T d given by 
Eq. OH]). The temperature profile is determined by the rate at which photons are generated, 
dn r /dt = riff + riDc, subject to n r = e rs /(3kT), with e rs given by the jump conditions. The 
equation governing the change in temperature across the thermalization layer is derived in 
appendix [X] (Eq. |A16j ). To obtain the temperature profile we numerically integrated Eq. 
(IA16P from T = T s at r d = up to T = T d at some distance r d = r dQ , after which the 
temperature remains constant at T d . Here r d is the distance from the shock front in units of 
the Thomson length. Examples are shown in figured] for upstream density rib u = 10 15 cm -3 
(corresponding to rj = 300 in (IA12[) and (IA13[) ) and different values of 7„. When the shock 
propagates in an expanding medium, as considered in §3, adiabatic cooling must be taken 
into account. The temperature then obeys Eq. ( 1A9j) . which is a generalization of f ]A16j) . 



2.2. NRRMS 

In the non-relativistic limit the jump conditions ([2])-(]4]) reduce to 

ft = ^(l + 8p/#), (12) 

18 

e rs = 3p rs = —n bu Plm p c 2 {\ -p/6#j), (13) 

where p = Pru/{nb u mpC 2 ) as before. It is readily seen that a shock can form provided 
> 4p/3. At lower velocities the upstream flow is subsonic. The advected photon density 
in the immediate downstream is given by n rs = (3 u n ru / (3 d = 7(1 + 8pf3~ 2 )~ 1 n ru . Substituting 
the latter result into ( ]T3l) we obtain the temperature in the immediate downstream: 

kT s = e rs /n rs ~ 10^4 + 8pf3~ 2 ) keV. (14) 

As noted above, the requirement that the flow is supersonic implies 8p/3~ 2 < 6, so that 
kT s < 7Qnl 1 f3 2 Ll keV. For the black body temperature we have 

kT d = 0.43n^ 5 /?y 2 keV. (15) 

The thermalization length is obtained as before. For free-free emission it is 

r'ff = 2 x lQ%} x n^^l\ (16) 

and for DC 

r' DC = 5.4 x lO^^n^/^l + 8p(3~ 2 ). (17) 
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3. Sub-photospheric shocks in GRBs 

We now use the above results to calculate the temperature downstream of a shock 
that breaks out of the photosphere of a GRB outflow. Consider a conical fireball having an 
isotropic equivalent luminosity L iso = 10 52 L 52 ergs s _1 . We assume that the fireball is ejected 
with an initial Lorentz factor T ~ 1 from a compact central engine of radius Rq = 10 6 i?6 
cm, and that it carries baryons with an isotropic mass loss rate Mj S0 . The properties of 
the fireball, and in particular the location of the photosphere, depend of the dimensionless 
enthalpy r\ = L iso / (M iso c 2 ). When r\ < r] c , where 

^(^r^ 8 * 1 ^ 1 ^ « i8) 

the fireball is sufficiently opaque, such that the radiation is trapped during the entire ac- 
celeration phase. The major fraction of the explosion energy is then converted into bulk 
kinetic energy of the baryons, and the fireball reaches a terminal Lorentz factor Too ~ 77 at 
some radius r coast ~ r]R Q /T , beyond which it continues to coast. The photosphere is located 
somewhere in the coasting region, at r p h > r coast . On the other hand, when 77 > r\ c the 
fireball will become transparent already during the acceleration phase, before reaching the 
coasting radius r = rj c Ro/T . The Lorentz factor in that case may be close to r] c (Nakar et 
al. 2005). 

We consider situations wherein rj < rj c . In that case the optical depth in the coasting 
region, at some radius r > r coast , can be expressed in terms of the fireball parameters as 



r r 



/ a T n b T- 1 dr = ^. (19) 

J r v r r 



The photospheric radius r p h can be found from the condition T(r p h) = 1, and it is readily 
seen that r ph = {t] c /rj) A r coast . Also, at the coasting radius r(r coast ) = (^ c /^) 4 . For the Lorentz 
factors inferred from observations, i] > 0.1r) c , the optical depth above the coasting radius, 
where shocks are likely to form, satisfies r 10 4 . 

The above results hold for a conical flow. Their generalization to a collimated flow is 
derived in appendix [Bj where it is found that for a flow geometry characterized by a cross- 
sectional radius a(z) oc z K , z ph /z coast = (fj c /vY 3+1 ^ K \ with fj c = T (r] c /T ) 4K ^ 1+3K \ For k < 1 
we have fj c < rj c , and the condition that the photosphere will be located in the coasting 
region occurs at lower values of 77. 

Now, suppose that the internal energy of the outflow at the coasting radius is dominated 
by radiation, as expected for hydrodynamic fireballs. Then, if the outflow remains adiabatic 
during the coasting phase, the internal energy changes as e r oc % , and the photon density 
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as n r oc rife. Consequently, the ratio n = n r /rib is conserved for an adiabatic flow. From 
f fT9|) one also obtains for the energy density of the radiation near the photosphere, e r {r ph ) = 
(i]/r] c ) 32 / 3 e T (r coast ), and for the flux r 2 h e r (r ph ) oc (v/r]c) 8/3 r 2 oast e r (r coast ). So unless 77 ~ r} c , 
or a considerable fraction of the bulk energy dissipates near the photosphere, the radiative 
efficiency of the fireball is extremely small. 

To estimate the photon density near the photosphere we suppose that the GRB outflow 
is adiabatic from its injection point at r = Rq up to the sub-photospheric region where shocks 
form. Since, as shown above, n is conserved along adiabatic streamlines, its value near the 
photosphere equals to its value at the injection point. From the hydrodynamic equations it 
can be readily shown that hT is conserved along streamlines, where h = 1 + 4p r / (n m p c 2 ) is 
the dimensionless enthalpy per baryon. Thus, h T = r], and since to a good approximation 
h = Ap r0 1 (n b0 m p c 2 ) we have@ 



3r] f m p c 



ft = i3Ti I ik ) ■ (20) 

The temperature T can be found using energy conservation: (IIolTq /3)rQ/3 c47ri?Q = L iso , 
and we adopt (3o = 1. Substituting T thereby obtained into Eq. (120]) gives 

h = 1.8 x 10 5 (77/77 c )^ /4 r o - 1/4 . (21) 

It is worth noting that the same result can be obtained upon assuming that at the coasting 
radius half of the outflow energy is carried by radiation. We emphasize that the dimensionless 
entropy given in Eq. (12"T1) depends only on the total power and baryon load of the fireball, 
and not its structure. It therefore holds for any outflow geometry. 

The ratio of photon pressure and rest mass energy density at a given optical depth r > 1 
can be obtained using ( !T9|) and (|20|) with T T Q = r]T coast : 

Pr( r r) Pr{r coast) ( r C oast \ 1 / n 2/3 /■_,/_, \ 8/3 



P(r) = P TJ 2 = P 7" 2 ™ - (TT /d (v/VcT /3 , (22) 
n h [r T )m v c l n b {r coast )m p c 2 \ r T J 

where r T = T~ x {r\ c jri)^r coas t is the radius at which the optical depth equals r. Thus, for 
i] << T] c , p << 1 near the photosphere. The corresponding temperature is readily found 
from Eq. (|23J) : kT(r) = 2.5r 1 / 6 (r//^c) 5/3 (r /i?6) 1/4 keV. 

Consider now a RRMS forming in the coasting region, at a radius r$ < r p h, where 
r p h is the photospheric radius, . As mentioned above, in its rest frame the shock has a 



2 At the base of the flow, where the temperature exceeds a few MeV, the radiation is in thermodynamic 
equilibrium with the e^- pairs, and the pressure is po = llaTg/12. 
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characteristic width A' s ~ (crn^) 1 , where n& u is the proper density of the unshocked gas. 
In the Lab frame the shock width is given by 

— = ^ = V-\r/r ph f = ^)-\ (23) 
r p h Wph 

where r = <JTn bu r/i] = (r p h/r). Evidently, the shock transition layer broadens as the shock 
propagates outwards. Since in the frame of the downstream fluid the shock velocity is 
(3d = 1/3, the comoving length of the shocked layer (i.e., the downstream region) at some 
radius r is A^ ~ r/(2r]), and its optical depth is Ar d = axribd A' d = V2 1u tE Consequently, 
the shocked layer becomes progressively more transparent as the shock evolves. The fraction 
of the dissipation energy contained inside the shock at a radius r is roughly Ar s / At^ ~ 
(a/27 u t) _1 , which can be significant near the photosphere. From Eq. (11 01) we have 

ATd ~ lO^C^VV/^^/To)- 1 ^, (24) 



T ff 



and from Eq. 09} 



Aid - 2 x 10- 5 (v/Vcf^ lu (R 6 /r )-^A DC1 . (25) 



T DC 



Consequently, thermalization can occur only for shocks that form at small radii, r/r p h ~ 10 4 . 
Since r coast = {r\ j r\ c )^r v ^ the latter condition requires r]/r] c < 0.1 for shocks that form in the 
coasting region, at r > r coast . 

The observed temperature of a shell that just crossed the shock is found from Eqs. (TTTT) 
and fl2Q: 

kT S:0b = r,kT s ~ 27 u /3 u L^ 4 r^ 2 J R- 1/2 MeV, (26) 

and it is seen that it is independent of the bulk Lorentz factor. This temperature may be 
different for different shells if 7„ changes during shock propagation. The corresponding black 
body temperature, defined as Td = (e rs /a) 1/,/l , can be expressed as 

kT d , ob = 7]kT d ~ 1.3 x 10-V^52 V4 (7^) 1/2 keV. (27) 

The maximum observed frequency behind the shock is hv max = 3kT 8i0 b ~ 67^ MeV for our 
fiducial model. 

The above considerations suggest that the evolution of a shock that propagates in an in- 
homogeneous, relativistically expanding medium is far more complex than that of an infinite 



3 Note that the optical depth r is defined in terms of the density of unshocked fluid, n bu , hence the factor 
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planar shock. Firstly, the shock velocity 7 U /3 U may not be constant during shock propagation. 
Secondly, as shown above, the shocked shell becomes more transparent as the shock moves 
outwards. This means that photons that were trapped in downstream shells that passed 
the shock at small radii, can re-interact with the shock at larger radii. As a consequence, 
the evolution of shocked shells may not be adiabatic. Thirdly, the shock width increases, 
suggesting that near the photosphere a significant fraction of the dissipation energy may be 
contained inside the shock, which may have important implications for the observed spec- 
trum. In fact, it is likely that some shocked shells that regain transparency will become 
part of the shock transition layer. Nonetheless, from fj24j) and f l25|) it is evident that photon 
generation occurs predominantly at small radii (large optical depths), on scales over which 
density variations and adiabatic cooling are small, so that the number of photons in the shell 
can be computed to a good approximation using the planary shock model. The increase in 
photon density by the thermalization process implies that even if the shocked shells are 
re-heated at larger radii, the average photon energy will be smaller than 3kT s given by Eq. 
( )26|) . as the energy is now shared by more photons. A complete treatment must account 
for the temporal evolution of the shock structure and the re-heating of the shocked shells 
in a self-consistent manner. Such an analysis is beyond the scope of this paper. Below, we 
present sample spectra computed under simplifying assumptions that will be outlined in the 
next section. 

Similar conclusions are obtained for NRRMS. From Eqs. fTl4|) . f JT5|) and f l2Tj) we find 



The thermalization depth is found to be smaller by a factor of about 10 to 100 than that of 
RRMS, depending on the value of (3 U . 



To illustrate some properties of the emitted spectrum, we shall (i) assume that the 
shock Lorentz ^ u is constant during the evolution, with the exception of one example (figure 
[5]); (ii) ignore possible reheating of the shocked shells near the photosphere; (iii) ignore 
the contribution of the shock transition layer (which might be important) to the observed 
spectrum. As explained above, the spectrum inside the shock is very hard, extending up to 
about , y u m e c 2 {rj^ u m e c 2 in the observer frame). So, our analysis below does not account for 
the portion of the spectrum above kT St0 f, s . 



kT S)0b = V kT s ~ l^L^ 4 ry 2 i?- 1/2 (l + 8p/3" 2 ) MeV. 
kT d:0b = V kT d ~ 1.6 x 10-yL 5 - 1/4 /3,y 2 keV. 



(28) 
(29) 



3.1. The spectrum emitted by a single shock 



Once the shock breaks out of the photosphere, an observer will start receiving radiation 
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from the shocked shells. Over time deeper and deeper shells will reach the photosphere 
and emit. Let t = be the Lab time at which the shock reaches the photosphere, and 
let x s denotes the distance of some shocked shell from the photosphere at t — 0. The 
radius at which that shell was formed (i.e., crossed the shock) is r s = r ph — 2rfx S) and the 
corresponding optical depth is t s . The shell will reach the photosphere at time t = x s /c, 
with a temperature T p h(r s ) given as a solution of Eq. (1A9j) . The radiation trapped in the 
shell has a Wien spectrum. At the photosphere, 

U t = ,,/c) = ( 

where e r (r ph ) = 2>yln bu (r ph )m p c 2 (r s /r ph ) a/3 = 2>yln bu (r ph )m p c 2 n a/J ' (see Eq. dA7j)), and we 
invoke a density profile n u oc r~ a , a > 1, for the unshocked flow. The isotropic equivalent 
spectral luminosity emitted from the photosphere is given by L^it, r p h) = 4ir 2 r 2 h I u (t = x s /c), 
and the time integrated spectral energy distribution (SED) by 

uK = f vL v {t,r ph )dt oc £ (j^^j r^'Me-^.W,. (31) 

Here To corresponds to the radius r at which the shock was initially formed. 

Figures HI display the time integrated SED, vE u , for a conical flow (a = 2), and 
different values of the remaining parameters, assuming 7„ = const. All curves are normalized 
to their peak values. As seen, the integrated emission exhibits a roughly thermal spectrum, 
with a peak energy hv peak ~ 3kT s ^ ob , where T s oh given by Eq. ( 126|) for RRMS and by Eq. (|28|) 
for NRRMS. The conceivable effect of reheating near the photosphere and the contribution of 
the shock transition layer itself are ignored in these examples. The latter is naively expected 
to lead to a relatively hard spectral component above the peak energy hv pea k = 3kT Sj0 b, up 
to about hv ~ rj^ u m e c 2 . 

The effect of adiabatic cooling is merely to change the slope of the spectrum below the 
peak, as seen in figure HI For a shock forming at an optical depth r >> 1, the time integrated 
SED at energies between hv pea k/TQ^ 3 and hv pea k is a power law with a slope 1 + 3/a (i.e., 
vE v oc z/ 1+3 / a )). The latter result can be derived using the following heuristic argument: 
The energy emitted from the photosphere scales as vL u (t,r p h) oc (r s /r p h) a ^ 3 , and the time 
integrated SED as vE v ~ r s \vL v {t,r ph )\ oc ri 1+a ^ 3 \ The peak energy of the spectrum 
emitted from the photosphere at time t = x s /c scales as hv = 3kT ph {r s ) oc (r s /r ph ) a ^ 3 . 
Thus, vE v oc v 1+3/a . 

The weak dependence of the emitted spectrum on the optical depth r is a consequence 
of the fact that the dominant contribution to the time integrated emission comes form 



u 3 e -hv/kT ph (rs) 



(30) 
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regions near the photosphere. The reason is that our assumption of constant shock velocity, 
JuPu =const, implies uniform dissipation, viz., dUdiss = rfe r (r s )dr s oc dm, where dm oc 
nb(r)r a dr = n^ir v h)r^ h dr is the mass of a shocked shell. Now, the optical depth across 
a shell of mass dm scales as dr oc n\,{r)dr oc r~ a dr. Thus, the mass enclosed below the 
photosphere satisfies m(r) oc 1 — T~ l '( a ~ l \ and it is seen that it is dominated by shells of 
moderate optical depth. Because the power emitted from each shell equals the dissipated 
power dUdiss modified by adiabatic cooling, it is clear the integrated emission is dominated 
by shells created just below the photosphere. This is the reason for the relatively small effect 
of adiabatic cooling seen in figure HI From the above it is also evident that at most a fraction 
{jth — of the dissipated energy can thermalize, where Ty, >> 1 is the thermalization 
depth. Consequently, thermalization is ineffective in case of a uniform dissipation profile. 
This statement is true in general, regardless of the specific dissipation mechanism. 

The caveat to the above argument is our assumption of a constant dissipation. A non- 
uniform dissipation profile might be established if the shock weakens as it propagates, as 
expected, e.g., when a fast shell collides with a much thicker, slow shell. In such a case the 
shape of the emitted spectrum can be significantly altered, and in particular the thermal 
peak is expected to be shifted to much lower energies. This is illustrated in figure EJ where 
the spectrum emitted from a shock moving at constant Lorentz factor, 7 n = 10, is compared 
with that emitted from a decelerating shock with 7 M (r s ) = 10(ro/r s ) 1//2 , where r$ = 10~ 2 r p h 
is the initial deceleration radius of the shock. The latter profile describes the dynamics of 
a blast wave propagating in a medium having a density profile n& oc r~ 2 , and is naively 
expected at times longer than the crossing time of the reverse shock. 

In all the examples depicted in figures EHU the portion of the spectrum below the peak 
is much harder than that of a typical Band spectrum. This reflects a generic shortage 
in production of soft photons by sub-photospheric shocks. The soft component may be 
produced subsequently by nonthermal processes at larger radii, after the shock has emerged 
from the photosphere and became collisionless, or by multiple shock emission (see below). 
On the other hand, the spectrum produced by a mildly relativistic shock (7^ £ 2) just below 
the photosphere, as those exhibited in figure [31 or by a decelerating relativistic shock that 
forms at a relatively large optical depth, as depicted in figure can nicely feet the thermal 
component observed in sources such as GRB090902B. We emphasize that the efficiency of 
mildly relativistic, sub-photospheric shocks can be quite significant. For instance, a shock 
propagating at 7„ = 2 can dissipate up to 30% of the bulk energy. 

A softer spectrum below the peak may be produced via multiple shock emission. The 
energy dissipated behind a shock moving at a velocity (3 U scales as fi\ in the nonrelativistic 
regime (Eq. |13]). and likewise the peak energy (Eq. [28]). Thus, for a uniform distribution 
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of shock velocities we expect vE v oc v. In the ultra-relativistic regime the same argument 
yields uE u oc u 2 . Complications of this simple situation can be envisaged, but it is generally 
expected that the integrated emission from multiple shocks with a range of strengths, or 
a range of fromation radii, might lead to a spectrum that mimics a Band spectrum. An 
elaborated treatment of multiple shock emission will be presented elsewhere. 

4. Conclusions 

We have analyzed the properties of subphotospheric shocks that form during the prompt 
phase of GRBs, and computed, under simplifying assumptions, the time integrated spectrum 
emitted as the shock emerges from the photosphere. We have shown that such shocks are 
mediated by Compton scattering of radiation advected into the shock by the upstream flow, 
in difference from shocks that propagate in a cold medium, e.g., shock breakout in super- 
novae, where photons are generated inside the shock transition layer by free free emission. 
We also argued that the scale of the shock, a few Thomson lengths, is vastly larger than any 
kinetic scale involved, so that particle acceleration by the Fermi process is highly unlikely 
in such shocks. The observed spectrum reflects the temperature profile downstream of the 
shock, with a possible contribution at the highest energies from the shock transition layer 
(BML11). 

Quite generally, in shocks that are dominated by photon advection the temperature 
just downstream of the shock transition layer, as measured in the shock frame, depends 
on the Lorentz factor of the upstream fluid and the ratio of photon and baryon densities 
advected into the shock. In GRBs, the latter ratio is determined by the fireball parameters, 
specifically, the outflow power and the baryon load. Since the immediate post shock region 
is photon starved, the temperature there is well in excess of the black body temperature. 
Photon generation by double Compton and free-free emissions in the downstream flow tends 
to thermalize the radiation. For typical GRB parameters, we find that the scale of the 
thermalization layer is about 10 4 - 10 5 Thomson depths, so that the radiation produced by 
sub-photospheric shocks that form in the coasting region is unlikely to reach full thermal 
equilibrium. 

The spectrum emitted by a single shocked shell that has reached the photosphere is 
a Wien spectrum, with a characteristic temperature roughly equals the immediate post 
shock temperature, modified by adiabatic cooling and thermalization. The time integrated 
spectrum is the sum over all shells, and reflects the temperature profile behind the shock. 
We have demonstrated that the emission from a single shock has a prominent thermal peak. 
Such episodes can naturally account for the thermal component observed in GRB090902B 
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and similar sources. The location of the peak depends on the velocity profile of the shock; for 
a uniform dissipation profile (i.e., 7„ = const) the thermal peak is located at hv pea k = 3fcT S)0 &, 
with T s>ob given by Eq. flU} for RRMS and by Eq. flU} for NRRMS. The peak energy is 
significantly reduced in case of a decelerating shock (figure E]). 

In general, the spectrum emitted by a single shock is very hard below the thermal peak, 
much harder than a typical Band spectrum. However, we qualitatively demonstrated that 
multiple shock emission can in principle give rise to a much softer component below the 
peak, with vF v oc z/ s , 1 < s < 2, for mildly relativistic shocks with a uniform distribution 
of velocities. Further work is needed to quantify the effect of multipole shock emission in 
more realistic situations, but the above reasoning suggests that multipole shock emission 
may mimic a Band spectrum. 

Unlike the case of an infinite planary shock, the structure of a shock that propagates in a 
relativistically expanding, nonuniform medium evolves with time. To be concrete, the shock 
width broadens, and the shocked fluid downstream becomes progressively more transparent 
as the shock moves outwards. This raises the possibility that the trapped radiation down- 
stream of the shock will re-interact with the shock as it approaches the photosphere, leading 
to a redistribution of the radiation energy. Moreover, we find that a significant fraction of 
the total dissipation energy is contained inside the shock as it breaks out of the photosphere. 
Consequently, the hard radiation produced inside the shock transition layer may significantly 
contribute to the high energy portion (above 3fcT s obs ) of the observed spectrum. A power 
law component extending up to about r]m e c 2 £ 100 MeV is conceivable. 
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A. Computing the temperature profile 

The change in the photon density n r downstream of the shock, owing to free-free and 
double Compton emissions, is governed by the equation 



d li (n r u' M ) =hff + h DC , 



(Al) 
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where is the 4-velocity of the shocked gas. Assuming n e = n p = n b d in the downstream 
flow, the specific photon generation rate by Maxwellian e-p Bremsstrahlung emission can be 
expressed as 

hff = a e a T cn 2 bd (kT/m e c 2 y 1/2 A f f, (A2) 

where A// = 10A//1 = [—Ei(hh / min /m e c 2 )]gff, Ei(x) is the exponential integral of x, and 
gff is the Gaunt factor. 

Under the conditions considered here the Compton y parameter exceeds unity for shocks 
that form at optical depth r £ 10. For y » 1 the cutoff energy hu min is determined by the 
requirement that every newly generated photon can reach the thermal peak by scattering off 
the hot electrons before getting re-absorbed. This yields Aff ~ 10 at temperatures kT d of a 
few keV (see, e.g., KBW09 and BML11 for details). The rate of double Compton emission 

iS 16 

h DC = — a e a T cn bd (kT/m e c 2 ) 2 n r A DC , (A3) 

7T 

here A^c = hi(/cT /hv m i n )gDc- For the ratio of the two rates we have, 

= ^(kT/m e c 2 fl 2 {n r /n bd ){A DC /A ff ). (A4) 
n ff 71 

Now, if the shock propagates in a stationary conical flow of a constant Lorentz factor, 
u° — cq » 1, then Eq. (IA1I) reduces to 

^^-{r 2 n r 7]) = h ff + h DC . (A5) 
r z ar 

Dividing the last equation by n& u , and recalling that n^r 2 = const and n bd = \f^^ u nh u -, one 
obtains 

— {n r /n bu ) = -8a el 2 u (kT/m e c 2 )- 1/2 Aff(l + h DC /h ff ) (A6) 

in terms of the optical depth dr = a T n bu dr /rj, where Eq. (IA2I) has been employed. It can 
be readily shown that the last equation holds not only for a conical flow, but essentially for 
any geometry, provided the opening angle is sufficiently small. Behind the shock the flow is 
adiabatic, so that the energy density of the radiation field satisfies e r oc n bd . If the density 
of the medium in which the shock propagates has a power law profile, n bu oc r~ a , then 
e r oc r~ 4Q / 3 assuming j u =const. For a shell that crossed the shock at a radius r = r s (t), 
where r s (t) is the shock trajectory, e rs is given by Eq. ()6]). At r > r s 

e r = e rs (r s /r) A ^ 3 = 2n bul 2 u m p c 2 {r s /r) a l\ (A7) 
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Now, on scales smaller than the thermalization length, the trapped radiation has a Wien 
spectrum with n r = e r /3kT. We can therefore write 

kT=^f-{r s /rTl\ (A8) 
3{n r /n bu ) 

Substituting Eq. ( 1A8I) into ( 1A6I) and using ( 1A4I) . we arrive at 

dT_oT = A RfWw Ts / T f/3 + (A9) 
dr 3r 

where T = T/T s , a — a/(a — 1) (a > 1), and 

A* = 12a e (m e /m p )A // (A;T s /m e c 2 ) 1 / 2 , (A10) 
k k = 2 x 10 3 7u (A;T s /m e c 2 ) 3/2 . (All) 

For our fiducial model, L 52 = T = R 6 = 1, Eq. f[2"6"j) yields (kT s /m e c 2 ) = 4(7^/77), and Eqs. 
(jATQjl and flATT]) reduce to 



A K = 1.7 x 10- 4 r ? - 1 / 2 7 - 1 / 2 A //1 , (A12) 
k r = 1.6 x lOS^r/' 3 / 2 . (A13) 

Let T ph {r s ) denotes the temperature of a shell that crossed the shock at some radius r s < r ph 
(r = t s > 1) as it reaches the photosphere (r = 1). Specifically, T ph {r s ) = T s T{r = 1), with 
T s given by Eq. f l2"B"j) . For a given r s , T ph (r s ) is obtained by integrating Eq. fl A9j) from r = r s 
to r = 1, subject to the boundary condition T(r = r s ) = 1. Equation (1A9j) holds as long as 
n r < aT 3 /k or equivalently aT 4 > e rs (r s /r) 4Q//3 . 

In the non-relativistic limit the temperature T obeys Eq. (IA9[) with replaced by 

^ = 3x lO-^-V^A^Cl + 8^- 2 )- 3/2 , (M4) 
2 x 10 4 /3 4 r/~ 3 / 2 (l + 8/5/3" 2 ) 3/2 . (A15) 



To compute the temperature profile behind an infinite planar RRMS, it is convenient 
to transform to the variable r d = o^n^x', where x' is the distance from the shock front, as 
measured in the frame of the shock. The temperature then obeys the equation 

(3 d ^- = -A R f 3 / 2 [l + k r T 3 % (A16) 



with A R and k r given by Eqs. flA12j) and (IA13I) . respectively. The boundary condition in 



this case is T(rd = 0) = 1. Likewise, for a NRRMS, T obeys the same equation with A R , k r 
replaced by A NR , k nr given above. 
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B. Upstream conditions for a collimated flow 

In §3] we determined the conditions upstream of a sub-photospheric shock assuming 
radial geometry for the unshocked flow. Here, we generalize the analysis to collimating flows. 
We consider a jet moving along the z-axis in a channel of varying cross-sectional radius a(z). 
For simplicity we adopt the profile a(z) = a (z/ z ) K . The jet is characterized by a total power 
Lj = nbhm p c 3 r 2 7ra 2 , and mass flux M = nbm p cT7ra 2 , where h is the dimensionless enthalpy 
per baryon. In the acceleration zone, where h » 1, the bulk Lorentz factor increases as 
r = To(a/ao). The transition to the coasting regime, where r = rj = Lj/Mc 2 , occurs at a 
distance z coast at which a coast = a(z coast ) = r]a /T : 



Zcoast 



^o(w/ao) 1/K = ^o(Wr ) 1/K , (Bl) 



and it is seen that the coasting radius of a collimated flow (ac < 1) is larger than that of a 
radial flow by a factor (r]/Top 1 ~ K '' K . 

The optical depth at a distance z from the origin can be expressed as 

/ \ n-i VrLjZ cr T L iso rfczQ 

r(z)^a T n b zT =^ = ^ = ^ (B2) 

where the isotropic equivalent power is defined as Li so = ALjz 2 /a 2 , and rj c is given by Eq. 
f|T8|) with R = zq. Thus, in the coasting region, where T = r/, the optical depth is given by 
Eq. fTTTJT) also for a collimating flow. The photosphere is located at z p h = z rj^/ {T rj 3 ). 

Let fj c denotes the critical load for which z coast = z p h- From the above result for z ph and 
Eq. (1B1I) we obtain 

77c = r (7/ c /r ) 4K/(1+3K) . (B3) 

For a radial flow (k = 1) Eq. (1B3[) reduces to fj c = rj c , as it should. For a collimating flow 
f\ c < i] c . The relative location of the photosphere is given by z p h/ z coast = (f]c/f]Y 3+1 ^ K \ thus 
the photosphere will be located in the coasting region provided r\ < fj c . 

The dimensionless entropy n, Eq. (I2U|) . depends only on the conditions in the vicinity 
of the central engine, not on the structure of the flow. Consequently, Equation (j2ip remains 
valid also in the general case of a collimated flow. However, the ratio of radiation pressure 
and rest mass energy density does depend on the structure of the flow. For the outflow 
profile adopted above we have 

P(r) ~ ~ (t) 2 / 3 (77/?7c) 2(1 ~ k)/3k - (B4) 

This difference merely affect the condition for formation of subrelativistic shocks, viz. /3 U > 
4p/3, not so much the structure of RRMS. 
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For illustration, we compare the above results with the numerical simulations reported 
in Lazzati et al. (2009). They analyzed the propagation of a GRB jet through the envelope 
of a 16 solar mass Wolf-Rayet progenitor star. For numerical convenience the jet is injected 
at a distance of R in j = 10 9 cm with a total power of Lj = 5.3 x 10 50 erg/s, opening angle 
#o = 10°, initial Lorentz factor Tq = 5, and dimensionless enthalpy ho = 80 for which 
7] = hoTo = 400. Fig 2 in their paper indicates that the transition to the coasting regime 
occurs at r coast ~ 10 11 cm, thus it seems that in this simulation the core jet is nearly 
radial above the injection radius Ri n j. For the above values we find an isotropic equivalent 
luminosity at the injection point of Li SO ~ 7 x 10 52 erg/s. It is difficult to predict what would 
the structure of the flow if injected at a realistic radius. If we naively interpolate the solution 
down to the central engine, at Ro = 10 6 Rq cm, with R$ ~ 1 for a 3 solar mass black hole, we 
obtain r] c = 3 x 10 3 Rq 1//4 . This yields a photospheric radius of r ph = (r)j/r) 3 )R ~ 1.2 x 10 12 
cm, independent of Ro, in agreement with their numerical result. For the channel profile 
adopted above, a coasting radius of r coast ~ 10 11 cm corresponds to a collimating jet with 
k = 0.6. Near the photosphere, where enhanced dissipation by collimation shocks takes place 
(see figure 2 in Lazzati et al. 2009), Eq. (1B4|) gives p ~ 0.1. Thus, NRRMS can form there 
provided (3 S > 0.4. The spectra emitted by such shocks are well represented by Fig 3. For 
RRMS our model provides a reasonable description even if form below the coasting radius, 
at larger optical depths. 
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Fig. 1. — Temperature profile in the thermalization layer of an infinite planar shock, for 
an upstream density = 10 15 cm -3 , and different values of the shock Lorentz factor 7 U 
(indicated by the numbers that label the curves). 
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Fig. 2. — Time integrated spectral energy distribution for rj = 200, r = 10, a = 2 (conical 
flow), and different values of the shock Lorentz factor 7„ (indicated by the numbers that 
label the curves). 
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Fig. 3. — Same as figure [2] for NRRMS with different shock velocities /3, 
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Fig. 4. — Same as figure |2] for r\ = 600, ^ u = 2, and different values of r , as indicated. The 
effect of adiabatic cooing is clearly seen at large values of r . 
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Fig. 5. — The solid line delineates the time integrated SED produced by a decelerating shock 
having a Lorentz factor profile 7 u (t s ) = \Q{t s /tq) 1 I 2 , with r = r ph /r = 100. The unshocked 
flow in this example is taken to be conical (a = 2) with 77 = 200. The spectrum produced by 
a shock moving at a constant Lorentz factor, 7„ = 10, is displayed for a comparison (dashed 
line). 



